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1.  Introduction 


Very  long  finned  projectiles  carrying  very  dense  metallic  rods  have  been  observed  to  be  subjected 
to  large  inelastic  defonnations  during  hypersonic  flight  (/).  Spark  shadowgraphs  of  these 
projectiles  have  shown  elastic  bending  motion  with  amplitudes  as  large  as  the  rod’s  radius  (2).  It 
has  been  conjectured  that  the  cause  of  the  inelastic  deformation  was  the  bending  loads  associated 
with  pitching  and  yawing  motion  at  resonant  spin  frequencies.  These  frequencies  would  be  near 
to  the  aerodynamic  frequency  of  a  rigid  projectile  or  the  elastic  frequencies  of  the  rod. 

Spinning  motion  at  a  resonant  frequency  can  be  caused  by  a  nonlinear  roll  moment  associated 
with  the  roll  orientation  of  the  fins  or  by  mass  asymmetry  due  to  damage  or  poor  construction  of 
the  projectile.  This  spin-yaw  lockin  mechanism  has  been  discussed  for  a  rigid  missile  by  several 
authors  (3-7). 

The  linear  flight  motion  of  an  elastic  missile  has  been  considered  in  references  (8-13).  In 
reference  (11),  a  very  simple  theoretical  model  of  a  projectile  composed  of  three  components 
connected  by  two  bent  mass-less  elastic  beams  was  used  to  approximate  an  elastic  missile.  For 
appropriate  selection  of  beam  parameters  spin-yaw  lockin  was  demonstrated  and  large  amplitude 
oscillations  occurred.  In  references  (12)  and  (13),  the  correct  linear  partial-differential  equation 
for  a  continuously  elastic  missile  was  derived,  and  special  solutions  for  harmonic  transient 
motion  and  trim  motion  were  obtained.  A  nonlinear  roll  moment  was  shown  to  be  capable  of 
producing  equilibrium  spin  near  resonance.  Although  the  lateral  motion  of  the  pitching  and 
yawing  missile  was  linearized,  it  was  necessary  to  retain  the  quadratic  terms  in  the  roll  equation. 
A  time  history  of  motion  going  to  lockin  was  not  obtained  at  that  time.  It  is  the  purpose  of  this 
report  to  compute  such  a  time  history  by  use  of  the  finite  element  method  (FEM)  (14). 

The  lateral  motion  of  a  symmetric  rigid  missile  has  often  been  described  by  complex  variables 
(15),  and  complex  variables  were  used  in  references  (8)  and  (11-13)  to  describe  the  lateral 
flexing  motion  of  the  rotationally  symmetric  rod.  The  equations  of  motion  of  a  spinning  missile 
are  usually  derived  by  vector  analysis  and  Newtonian  mechanics,  but  FEM  equations  for  elastic 
bodies  require  the  use  of  the  more  sophisticated  Lagrangian  mechanics.  The  Lagrangian  of  a 
system  must  be  stated  in  body-fixed  coordinates.  Separate  differential  equations  in  the  two 
lateral  directions  can  be  obtained  from  the  appropriate  Lagrangian.  Pairs  of  equations  are 
combined  and  complex  variables  introduced  to  yield  half  as  many  complex  equations. 

If  five  elements  are  used  to  describe  the  pitching  and  flexing  motion  of  the  missile,  20 
parameters,  usually  called  connectors,  are  introduced  and  25  second-order  differential  equations 
and  one  first  order  differential  equation  are  required.  The  introduction  of  complex  angles  and 
complex  connections  reduce  this  system  to  12  second-order  complex  differential  equations  and 
two  real  differential  equations. 
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2.  Coordinate  System 


The  elastic  missile  is  assumed  to  consist  of  a  very  heavy  elastic  circular  rod  of  fineness  ratio,  L, 
and  mass,  m,  embedded  in  a  very  light  symmetric  aerodynamic  structure  that  may  be  longer  than 
the  rod.  The  rod’s  axial  moment  of  inertia  is  Ix  and  its  transverse  moment  of  inertia  about  its 
center  is  Ito.  The  rod’s  diameter  can  vary  over  its  length,  and  its  maximum  diameter  will  be 
denoted  by  “d.”  All  distances  will  be  expressed  as  multiples  of  the  rod  diameter  and  its  length  is 
Ld.  A  nose  windshield  of  length  x23d  may  be  attached  to  the  forward  end  of  the  rod,  and  the 
fins  may  extend  beyond  the  end  of  the  rod  at  a  distance  of  x01d  .  Thus,  the  rod  is  located 
between  xt  =  -L/2  and  x2  =  L/2  ,  while  the  aerodynamic  structure  extends  from  x0  =  x  -x01 
to  X3  =  x2  +  x23 . 

An  earth-fixed  coordinate  system  will  be  used  with  the  Xe-axis  oriented  along  the  initial 
direction  of  the  missile’s  velocity  vector.  The  Ze-axis  is  downward-pointing  and  the  Ye-axis  is 
determined  by  the  right  hand  rule.  A  nonrotating  coordinate  system,  XYZ  is  then  defined  with 
origin  always  at  the  center  of  the  rod  and  the  X-axis  tangent  there.  The  X-axis  pitches  through 
the  angle,  0  ,  and  yaws  through  the  angle,  \|/ ,  with  respect  to  the  Xe-axis.  Body-fixed  coordinates 
XYb-Zb  are  now  defined  for  which  the  Yb-Zb  axes  rotate  with  the  missile  through  the  roll  angle, 
cj) ,  with  respect  to  the  Y-Z  axes. 

We  will  conceptually  slice  the  missile  into  a  large  number  of  thin  disks  perpendicular  to  the 
X-axis  with  thickness,  dx.  When  the  rod  flexes,  the  disks  shift  laterally  perpendicular  to  the 
X-axis  and  cant  to  be  perpendicular  to  the  centerline  of  the  disks.  This  canting  action  neglects 
the  shear  deformation  of  the  rod,  and  this  constraint  is  called  the  Bernoulli  assumption  (14).  The 
lateral  displacement  of  a  disk  has  body-fixed  coordinates  8by  and  5bz ,  and  the  disk  is  canted  at 

angles  Tby  and  Tbz. 


r 


by 


fix 


=  55^_ 

fix 


(1) 


It  important  to  note  that  at  the  central  disk 

s„  (o,t)= sb,  (o.t) = rby  (o,t)=rbz(o,t)=o .  (2) 

Reference  (12)  used  the  nonspinning  elastic  coordinate  system  with  XYZ  axes  (c()  =  0) .  The 

lateral  displacements  of  a  disk  in  this  elastic  coordinate  system  are  shown  in  figures  1  and  2  and 
can  be  computed  from  body-fixed  quantities. 

SE=SE,+iSEl=(Sby+iSbz)  e'*,  (3) 


and 


2 


Figure  1.  X-Z  coordinates  of  cross-sectional  disk. 


Figure  2.  X-Y  coordinates  of  cross-sectional  disk. 
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r  =  r,+irl=(rby  +  irte)  e». 


(4) 


The  earth-fixed  coordinates  of  the  central  disk  are  xe,ye,  and  ze ,  and  the  earth-fixed  coordinates 
of  the  other  disks  are  computed  in  terms  of  the  central  disk  earth-fixed  coordinates,  their 
body-fixed  displacements,  and  the  Euler  angles  9,  \j/,  and  (j) . 


and 


x ,  =  x 

de  < 


+  x 


l-(v2+02)/2]-vRe{(5by  +  i5bz)e^}  +  0Im{(5by+i5bz)« 


yde  =  ye  +  XV  +  Re  j(sby  +  i8bz )  e^ } , 

zde  =ze-x0  +  Im{(8by+i6bz)e4}. 


(5) 

(6) 

(V) 


ye,  ze,  0,  \| /,  8by,  8bz  and  their  derivatives  are  assumed  to  be  small  quantities,  but  xe  and  xe  are 

not  small.  Lagrangian  dynamics  yields  the  correct  linearized  differential  equations  for  these  six 
variables  when  the  kinetic  energy  expansion  retains  all  quadratic  terms  in  these  variables.  Thus, 
equations  6  and  7  consist  of  only  linear  terms,  but  equation  5,  which  contains  xe,  retains 
quadratic  terms  in  these  variables.  The  coefficient  of  x  in  equation  5,  for  example,  is  the 
quadratic  form  of  the  cosine  of  the  angle  between  the  Xe-axis  and  the  X-axis. 

The  angular  velocity  of  the  central  disk  in  the  body-fixed  coordinates  is 


P  =  (j)-Xj /0, 

(8) 

qb  =Re{(0  +  iijr)  e^} , 

(9) 

and 

rb  =  fin  j(0  +  i\j/)  e-1*  j  . 

(10) 

The  angular  velocity  of  any  other  disk  along  axes  aligned  with  the  disk’s  axes  of  symmetry 

is 

Pd  =  p  (1  -  rr/2)  +  qarby  +  rarte  =  P  +  Re }(e + hj<  -  *r/2 + if )  r} , 

(ID 

q<  =  qb  -fte  —  4>rby  =  Re{(e+hj<+if )  e-*} , 

(12) 

and 

r,  =  q.  +  r.,;.  ■  =  im  j(o+iij/  +  ir)  iT4}. 

(13) 
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3.  Motion  of  Central  Disk 


The  mass  of  each  circular  disk  is  (m/L)  p,  (x)  dx ;  its  roll  moment  of  inertia  is 
2ad  (md2)  p2(x)  dx ,  and  its  transverse  moment  of  inertia  is  ad  ( md2)  p2(x)  dx  ,  where  ad  is 
( 1 6L)  1 .  pj  (x)  and  p,  (x)  describe  the  variation  of  mass  and  moments  of  inertia  along  the  rod, 
and  both  become  unity  for  a  homogeneous  rod  with  constant  diameter.  The  kinetic  energy  of  a 
disk  is  therefore 

Tddx  =  [(md2Pl  /2L) (x2e  +  y2e  +  z2e )  +  (admd2p2 /l) (2p2  +  q2  +  r2 )]  dx  •  (14 

The  total  kinetic  energy  of  the  missile  can  be  obtained  by  integrating  Td  over  the  length  of  the 


T  =  J  Tadx  =  T„  +  md2  (T„  +  T„ )  +  (md2/2L)  (T;|  +  T,, ) , 

X1 

where  TR  =(md2/2)(xe2 +  y2 +22)  +  Ixp2/2  +  It0(\j/2 +02)/2 
+  md2xc  [ye\j /  -  ze0  -  xe  (w  +  00)] , 

Tn  =-xeRe{(\|/  +  i0)8c+(¥  +  i0)8c}5 
Ti2  =Re{(ye-iZe)5c+(\i/  +  i0)(j6+J8)}, 


T21  =  j  (S2y  +52Ez)p1dx  +  (adL)  J  (f ;  +  f  2)p2dx, 


T22=2cj,(adL)Re  J(2if-cj,r)rp2dx 


A 2 

xc  =(VL)  Jxpidx. 


The  ye  and  ze  components  of  the  central  disk  velocity  can  be  approximated  by  linear  relations  in 
angles  of  pitch  and  yaw  with  respect  to  inertia  axes  (0,  \(/)  and  angles  of  attack  and  sideslip  with 
respect  to  the  velocity  vector  (a,P)  as  shown  in  figures  1  and  2.  The  magnitude  of  the  velocity 
vector  is  V. 


ye  =(V/d)  (P+y): 
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Ze=(V/d)  (a-0) . 


(17) 


Equations  16  and  17  can  be  written  as  a  single  complex  equation: 

ye+iZe=(V/d)  (qj+qJ,  (18) 


where  =  P  +  ia  and  qle  =  \|/  -  i0 . 


In  reference  (12),  the  linear  aerodynamic  force  loading  is  expressed  in  terms  of  three  force 
distribution  functions,  cD  (x),  cfl  (x),  and  cf2  (x)  and  the  base  pressure  coefficient,  Cobp,  plus  a 

body-fixed  force  associated  with  possible  bent  fins.  Because  the  lateral  motion  of  the  missile  is 
quite  small,  q1  =  -qle  and  the  aerodynamic  damping  force  terms  in  the  aerodynamic  loading  on 

the  aerodynamic  structure  can  be  combined.  This  aerodynamic  loading  in  nonrotating  elastic 
coordinates  is 


^p  =  -g,cD(x), 

dx 


dfy  .  dF 

— -  +  i — - 

dx  dx 


•'n 


(x)  [q1-r-rBF(x)ei*+(8E-xq1)  (d/V) 

+  cf2(x)  (2q, -f -i<j>rBF(x)  e'*)  (d/V) 


and 


^xbp  Sl^Dph  • 


(19) 

(20) 


(21) 


The  total  aerodynamic  force  acting  on  the  aerodynamic  structure  is  given  by  the  integrals  of 
equations  19  and  20  and  by  adding  the  base  drag  of  equation  21  to  the  axial  force, 


Fx  =  -&CD 


(22) 


and 

F  =  -g,[c1qi+cj('i,‘1/V)-CNBFei‘-J1(t)-i2 (t)  (d/V)],  (23) 

where  various  functions  are  defined  in  appendices  A  and  B. 

The  primary  components  of  drag  are  head  drag  and  base  pressure  drag.  The  third  component  is 
skin  friction  drag  that  is  ~  15%  of  the  total  drag  and  will  be  neglected  in  this  report  to  simplify 
the  FEM  calculations. 

Similarly,  the  total  aerodynamic  moment  about  the  rod’s  center  can  be  computed  from  the 
transverse  aerodynamic  force  and  a  small  axial  force  contribution: 
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(24) 


M  =  My+iMz 

=  -i(g1d)[c3q1+c4(qId/V)-CMHFei*-J3(t)-j4(t)(d/V)-J5(t)]. 

The  motion  of  the  central  disk  is  described  by  the  variables  xe,  ye,  ze,  (j),  \|/,  and  9 .  The  motion 
of  any  other  section  on  the  aerodynamic  structure  located  at  x  with  thickness  dx  is  the  sum  of  this 
motion  plus  its  motion  relative  to  the  central  disk.  The  work  done  on  any  section  caused  by 
motion  of  the  central  disk  is 

(dWcd)x  dx  =  dQIxdxe  +dQlydye  +dQlzdze 

+  dQ2xd(j)  +  dQ2ydqley  +  dQ2zdqlez ,  (25) 

where  the  sectional  generalized  forces  are  defined  in  appendix  B. 

The  total  generalized  forces  can  be  obtained  by  integrating  the  sectional  generalized  forces  and 
adding  the  work  done  by  the  base  pressure  drag.  These  are  also  given  in  appendix  C. 

The  linearized  Lagrangian  differential  equation  for  xeand  (j)  gives  the  usual  drag  and  spin 
equation, 

mV  =  mxe  =  -gjCD ,  (26) 

and 

M>  =  MX.  (27) 

The  Lagrangian  differential  equation  for  ze  can  be  multiplied  by  i  and  added  to  the  Lagrangian 
differential  equation  for  ye  to  yield  a  single  second-order  differential  equation  in  complex 
variables.  Next,  the  Lagrangian  differential  equation  for  0  can  be  multiplied  by  i  and  subtracted 
from  the  Lagrangian  differential  equation  for  \j / , 

m  [V  (q,  +  q2 )  +  Vq,  +  8cd  +  q2xcd]  =  F  (28) 

and 

(lt0+md2x2)  q2-i(j)Ixq2+md2[j6+j8  +  xc5c]  =  -iM-xcFd,  (29) 

x2  x2 

where  xc  =(1/L)  J  xptdx  ,  5C  =(l/L)  J  S^dx  ,  and  q2  =  qle. 

X,  Xj 

These  differential  equations  for  xc  =  0  are  the  same  as  those  derived  by  Newtonian  mechanics  in 
reference  (12).  Equations  28  and  29  contain  eight  integrals  of  SE  and  T  .  Three  of  these  are 
integrals  of  dynamics  properties  (5C,J6,J8)  over  the  rod  (x,,x2)  and  five  are  integrals  over 
aerodynamic  properties  ( J, ,  J 2  ,J3  ,J4,J5)  over  the  aerodynamic  structure  (x0,x3). 
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For  constant  spin  and  velocity,  a  rigid  unbent  missile  with  its  center  of  mass  at  the  rod  center, 
equations  28  and  29  predict  a  simple  epicyclic  angular  motion: 

q,  =  K10eM+i<”  +K20eM+4\  (30) 

where  <j)m  =  cj>Ix /2It  ± ^-(g,d/It)  c3  +((j)Ix/2It) 2  . 

For  zero  spin,  <j>j  =  -<j)2  =  coR  ,  where  coR  =  J(g1d/It)|c3|  . 

4.  FEM 

The  rod  is  assumed  to  be  represented  by  the  sum  of  an  inelastic  bent  component  rotating  with  the 
missile  and  an  elastic  deformation, 

5E(x,t)  =  SEB(x)  e4  +  5E(x,t)  ;  Xj<x<x2,  (31) 

and 

5b(x5t)  =  5EB(x)  +  8b(x,t)  ;  Xl<x<x2,  (32) 

,  x  dSFR (0) 

where  8FB  ( 0)  = - - —  =  0  . 

EB  w  dx 

Because  the  aerodynamic  structure  is  rigidly  attached  to  the  rod, 

5E(x,t)  =  8E(x1,t)  +  (x-x1)  r(xj,t)  ;  x0<x<x!,  (33) 

and 

8E (x,t)  =  8E (x2,t)  +  (x-x2)  r(x2,t)  ;  x2<x<x3.  (34) 

The  motion  of  the  elastic  component  of  the  rod  is  controlled  by  the  elasticity  of  the  rod  and  the 
aerodynamic  force  acting  on  it. 

FEM  is  a  very  powerful  method  for  calculating  the  time  history  of  the  elastic  flexing  motion. 

The  rod  is  divided  into  nj  elements  of  length  Le  =  Lj n  .  The  shape  of  the  j-th  element  is  given 

by  a  linear  combination  of  third-order  Hermitian  polynomials  (see  appendix  C). 

5by(x,t)  =  X9bPy(t)  Np(Z)’  (35) 

i 

and 
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(36) 


Ux,t)  =  XqbPz(t)Np(z), 

1 

where  x  =  Le (z  +  Zj);  zj=x,/Le+j-l  0<z<l. 

The  coefficients  of  the  polynomials  are  functions  of  time  and  are  called  connectors.  The  first 
two  are  the  deflection  and  slope  of  the  left  end  of  the  element  and  the  third  and  fourth  are  the 
deflection  and  slope  of  the  right  end.  To  ensure  continuity  in  deflection  and  slope  at  junction 
points,  the  corresponding  pairs  of  connectors  are  equal.  We  will  consider  only  an  odd  number  of 
elements  and  require  that  the  connecters  of  the  central  element  satisfy  equation  2. 

4b3  —  5qbl  +  Leqb2 ,  (3V) 

and 

qb4  =  24L;1qbi  +  5qb2>  (38) 

where  qbq  =  qbqy  +  iqbqz  • 

For  n  elements,  there  are  2ir  independent  complex  connectors.  It  is  convenient  to  let  the  index 

for  the  connectors  run  from  3  to  nt,  where  nt  =  2nj  +  2.  The  usual  FEM  procedure  calculates  the 
elastic  parts  of  the  integrals  in  equations  28  and  29  by  first  obtaining  integrals  over  each  element. 
These  are  linear  functions  of  that  element’s  four  connector  functions.  Due  to  equations  37  and 
38,  the  central  element  has  only  two  independent  connector  functions,  and  the  next  adjacent 
element  connector  functions  are  related  to  them.  Integrals  for  these  elements  are  specially 
computed.  The  desired  elastic  integrals  are  sums  of  these  sub  integrals  and  are  linear 
combinations  of  the  2n.  complex  connector  functions.  The  three  integrals  can  be  written  as 

St=L-'f]a,11qn+8CBe1*,  (39) 

3 

J,=b'Za!.q.+J»e",  (40) 


and 


J8  =adZb2„(qn-24q„)-4J8Bi^,  (41) 

3 

where  qn  =qhnc1'1'. 

The  five  integrals  have  contributions  from  the  aerodynamic  structure  extensions  as  well  as  from 
each  element, 
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(42) 


J,  +  h  (d/V)  =  X[(f,„  +  f,„  )  q„  +(g,„  +  g.,„)(q„d/V)] 

3 

+  [j,B  +  i(4.d/v)j;B]e'*, 

J,  +  i4  (d/V)  -  £[(4,  +  f,Jn )  qn  +(g!n  +  g,2„ )  (q„d/V)] 

3 

+  [j3B+i(<t>d/v)j4B]ei», 
and 

J5=Z(-CDa,„L-'+h„)qn+J5Bei).  (44) 

3 

Expressions  for  ScB,  JkB,  han  and  all  FEM  coefficients  in  equations  39-44  for  1,  3,  5,  and  7 
elements  are  given  in  appendices  A,  D,  and  E. 

Equations  39-44  can  be  used  with  equations  28  and  29  to  write  these  two  complex  differential 
equations  in  a  standard  fonnat: 

Z  [Rm„qn  +  (smn  +  ifc  )  qn  +  (Tm  +  i(j)T*n )  q„  ]  =  tm  exp  (ic|>)  .  (45) 

n=l 

The  coefficients  for  m  =  1,2  are  given  in  appendix  F. 

5.  Flexing  Motion 

In  order  to  derive  the  equations  for  flexing  motion,  the  kinetic  energy  given  by  equation  1 5  must 
be  expressed  in  tenns  of  the  rigid-body  parameters,  plus  the  connector  functions.  The  usual 
FEM  process  is  used  to  express  T2j  in  tenns  of  two  matrices,  two  vectors,  and  two  constants. 

T21  =ZZkmn[  qmyqny  +  qmzqzn]-2Re(4e^ZkBmqm]  +  (j)2IxB1L,  (46) 

3  3  [  3  J 

and 

T22  =  2<j»dL  Re  jz  Z  bmn  (2iqm  -  ‘km  )  4n  +  2hT'*Z  ^ bBm  (qm  +  2ikm  )  j  +  ¥Kb2L  »  (42) 

where  IxHI,IxB2  are  defined  in  appendix  A  and  kmn,klin,binn,bHn  are  defined  in  Appendix  G. 

The  potential  energy  stored  by  an  elastic  deformation  is 
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P-E-  =  (1/2)  |  (EI/d) 

X1 

where  E(x)  l(x)  =  E0I0p3 (x) . 

E0  and  I0  are  values  of  Young’s  modulus  and  the  area  moment  of  inertia  at  the  rod’s  center. 
p3  (x)  gives  the  variation  of  their  product  along  the  rod.  For  a  homogeneous  rod  with  constant 

diameter,  the  product  is  Eolo.  The  potential  energy  integral  can  be  replaced  by  a  sum  of  integrals 
over  each  of  the  m  elements.  After  integrating  over  each  element,  the  results  can  be  summed  to 

yield  a  linear  combination  of  the  4nj  real  connectors, 

nt  ni 

V=(  md2COo  /2L)  XXKmyqbny+qbmzqbnz)  Cmn  ,  (49) 

m=3  n=3 

where  cOq  =  E0I0/md3 ,  cmn  are  defined  in  appendix  G. 

The  work  done  on  the  rod  by  this  elastic  force  can  be  expressed  in  terms  of  generalized  force 
terms, 

dWE  =  (  md2/L)  <*±[  ^Emy^^lbrny  Q  Emz  dq 

bmz  J  ?  (50) 

m=3 

nt 

where  QEmy  +iQEmz  =  ~YjCmnqhn  . 

n=3 

The  Kelvin-Voight  elastic  damping  force  (16),  which  is  proportional  to  the  time  derivative  of  the 
elastic  shear  force,  has  generalized  force  terms: 

dWD  =  (  md2/L)m;^[ 

QDmy^qbmy  Q  Dmz  dq 

bmz  J  ?  (51) 

m=3 

nt 

where  QDmy  +iQDmz  =  -[  2^)  £cmAn  >  ®i  =  (4.730/L)2  co0  . 

n=3 

The  scale  factor  2of 1  is  selected  so  that  k  =  1  corresponds  to  critical  damping  of  the  first  elastic 
mode. 

The  external  aerodynamic  force  is  divided  between  the  force  acting  on  the  rod  and  the  force 
acting  on  the  aerodynamic  structure  extensions  at  the  ends  of  the  rod.  The  total  work  done  on 
the  missile  is  the  sum  of  the  work  done  by  these  forces. 


<9  8 


by 


<9x2 

v  j 


+ 


<9  8 


bz 


<9x 


dx , 


(48) 
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(52) 


dWa=£[dWAj+dW„i' 

j  =  l 

nt 

=(gid)Z[  Q  Amy  bmy  Q  Amz  dq 

bmz  J? 

m=3 

nt 

where  QAmy  +  iQ^  =  X [(f™  +  L )  fin  +  (gmn  +  gamn ) fin  ]  e'1*  +  fBm  +  faBm  +  4 ( gBm  +  gaBm ) 

n=3 

and  L.gm. fa,n„ 5 gamn 5 fBm , gBm , faBm , gaBm are  defined  in  appendix  G. 

Equations  50-52,  in  conjunction  with  the  kinetic  energy  defined  by  equations  15,  46,  and  47,  can 
be  used  to  derive  the  4n  Lagrangian  differential  equations  for  the  q bmy ,  q hm/  ’s  for 

m  =  3,  4, ...2^.  +2  .  The  equation  for  qbmz  should  be  multiplied  by  i  and  added  to  the 
corresponding  equation  for  qbmy  to  yield  2m  complex  differential  equations  for  the  complex 
variables  qbm  .  If  the  qbm  are  replaced  by  qm  in  each  equation,  the  result  has  the  standard  form 
given  by  equation  45,  where  the  14nj  complex  coefficients  for  m  =  3,  4...2nj  +2  are  given  in 
appendix  F. 


6.  Special  Solutions 


A  rigid  symmetric  finned  missile  flying  with  constant  spin  has  two  natural  frequencies:  <j)1R  and 
<j)2R ,  where  <j)1R  =  -<j)2R  .  For  zero  spin,  each  of  the  flexure  frequencies  would  give  rise  to  two 
coning  frequencies,  ±coK .  The  frequencies  present  in  the  motion  of  an  elastic  projectile  would 
form  an  infinite  sequence,  where  the  first  two  would  be  related  to  <j)1R  and  <j)2R ,  while  the  later 
ones  would  evolve  from  ±coK,  i.e.,  <j)2K+1  =  coK  ;<j)2K+2  =  -coK  .  The  odd  numbered  modes  have 

positive  frequencies  and  are  called  positive  modes,  while  the  even  numbered  modes  are  called 
negative  modes.  For  nonzero  spin,  the  pairs  of  frequencies  would  no  longer  be  equal  in 
magnitude. 

Transient  harmonic  solutions  of  the  homogeneous  constant  spin  version  of  equation  45  have  the 
form 

q„  =  q-e"  A  -  Ak  k  =  1,  2,  3, ...2  (2nj  + 1) .  (53) 

After  substituting  equation  53  into  equation  45,  a  set  of  linear  homogeneous  algebraic  equations 
is  derived,  which  is  specified  by  an  nt  x  nt  matrix, 
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(54) 


u™=A!Rm+A  (S„+i+S'„)  +Tm+iiTl. 

For  nj  elements,  there  are  2nj  +  1  pairs  of  positive  and  negative  frequencies.  The  lowest 
frequency  pair  is  related  to  the  two  aerodynamic  frequencies  and  the  other  2nj  pairs  are 
approximations  of  first  2n,  elastic  frequencies  of  the  infinite  set  of  frequencies  for  the  elastic  rod. 

The  paired  frequencies  are  equal  in  magnitude  when  spin  is  zero  but  the  positive  frequency  has  a 
larger  magnitude  when  spin  is  positive.  For  an  elastic  beam  considered  by  Geradin  and  Rixen 
(14),  the  first  nj  approximations  are  close  to  the  correct  value.  This  accuracy  is  determined  by 
the  ability  of  nj  elements  to  describe  the  corresponding  mode  shape  correctly.  Geradin  and 
Rixen  (14)  also  observe  that  the  approximations  are  upper  bounds  of  correct  value. 

For  constant  spin,  the  trim  motion  in  response  to  the  spinning  bent  missile  forces  can  be  obtained 
by  assuming  a  solution  of  the  form 

q„  =  qnTe*  ■  (55) 

The  response  parameters,  qnT ,  vary  with  particular  values  of  the  constant  spin  and  are  solutions 
of  a  set  of  nt  linear  inhomogeneous  algebraic  equations, 

nt 

XBmnq„T  =tm  >  (56) 

n=l 

where  B  =  — cb2  (  R  +S*  )  +  i<j>  (  S  +T*  )  +  T  . 

In  reference  (12),  both  these  special  solutions  were  computed  from  iterative  solutions  of  a 
second-order  differential  equation  specified  by  special  boundary  conditions.  The  FEM  approach 
obtains  the  same  results  and  more  by  much  simpler  matrix  operations. 


7.  Quadratic  Roll  Equation 


The  linear  roll  equation  (equation  26)  was  derived  from  a  kinetic  energy  function,  which 
neglected  cubic  terms  in qly,  qlz,  q2y,  q2z,  qb  ,  andqbnz-  In  order  to  derive  the  quadratic  roll 

equation,  cubic  terms  involving  (j)  and  (j)  must  be  retained.  The  resulting  roll  equation  is  the  same 
as  that  derived  in  reference  (12)  from  Newtonian  mechanics, 

+  Re {md2  ( j7  -  iq2  J6  +  iq2  Js )  - iF6cd}  =  Mx ,  (57) 

where  Jy(t)  is  defined  in  appendix  F. 
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The  aerodynamic  roll  moment  is  the  X-component  of  the  aerodynamic  moment  about  the  center 
of  the  central  disk.  The  linear  roll  moment  coefficient  for  a  rigid  finned  projectile  is  usually 
expressed  in  terms  of  a  roll-damping  coefficient  and  a  steady  state  spin, 

(c<U=c-p(<i>-p.)  (d/vf  <58> 

The  steady  state  spin  is  usually  determined  by  a  differential  canting  of  the  fins  caused  either 
intentionally  by  the  designer  or  unintentionally  by  damage  to  the  projectile. 

The  roll  moment  of  one  of  the  projectile’s  thin  disks  is  the  sum  of  the  linear  roll  moment  it  has  as 
part  of  a  projectile  and  the  quadratic  roll  moment  induced  by  the  transverse  aerodynamic  force 
acting  on  its  lateral  displacement  relative  to  the  center  of  the  central  disk.  If  we  retain  only  the 
dominant  Cfi  term  in  equation  20,  the  quadratic  roll  moment  has  the  fonn 

(dMpW=[(<lFI)SE,-(dFy)  5Ee] 

=  (g,d)  cfl  Re  {  i[(q, -r-rBFe‘*)  +  (5E-xq,)  (d/V)]  SEjdx-  <59) 

The  total  aerodynamic  roll  moment  acting  on  the  projectile,  therefore, 

M.  =(g,d)  [  (C,),_+Re{  Qa(t)-ig-'F5c}  ],  (60) 

where  QA(t)  =  ijcn[(q1-r-rBFe,‘)  +  (5E-xq1)  (d/V)]  (o,  -S  )  dx . 

x0 

The  nonlinear  roll  moment  from  equation  60  can  be  placed  in  the  spin  equation  57  to  yield: 

I.x4  +  Re {md2  (Qd  - iq2  J6  +  iq2  J8 )  - g,d QA }  =  gld (Cf  )linear ,  (61) 

where  the  quadratic  terms  IX.QD  are  defined  in  appendix  F. 

For  trim  motion,  (j)  =  Re  {  QD  }  =  0,  q2  =  icj)q2Te“,>t  and  equation  61  become  a  simple  equality  of 
two  functions  of  p. 

f*  (♦)=*;(♦)•  (62) 

where  f,  =<j>-p„,  f,  =  -R{q„  -  (q2Tmd/g, )  (<j>76T  +iJST)}(Clpd/v)  .and 

x3 

Qat— ijcfi  (diT  —  —  rBF )  +  i^8ET  —  xq1T ) ^(j)d/vj  8Tdx. 

x0 

8et,Tt  are  computed  from  the  solution  to  equation  56  and  J6t,  -fix  are  J6,  Jg  evaluated  for 
8ET,rT .  For  moderate  spin  Qat  is  the  dominant  part  of  f).  Equilibrium  values  of  spin  are 
determined  by  the  intersection  of  these  two  curves.  Lockin  occurs  at  stable  equilibrium  spin. 
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Spin  lockin  can  only  occur  when  the  missile  has  some  rigid  asymmetry  such  as  that  described  by 
8eb  or  rBF  .  For  a  rigid  missile  SET  and  Tt  are  SEB  and  TB .  q1T  is  a  function  of  spin  which  has 

a  large  resonant  amplitude  when  spin  is  near  coR  .  For  an  elastic  missile  q1T,8ET,rT  are  all 
functions  of  spin  and  have  large  resonant  amplitudes  when  spin  is  near  <j)2k_l  or  near  <j>,k  for 
negative  spin. 


8.  Numerical  Results 


In  reference  (12),  calculations  were  made  for  a  20-cal.  fin-stabilized  rod,  flying  at  6000  ft/s.  The 
finned  missile  has  a  1-cal.  nose  extension  and  a  1-cal.  fin  extension  (figure  3).  The  mass  and 
aerodynamic  parameters  for  the  finned  missile  are  given  in  appendix  H.  The  bent  rod  will  be 
described  by  a  pair  of  quartic  curves. 


SEB  =  dux2  +d21x4  -10<x<0 

=  d]2x2  +  d22x4  0  <  x  <  10.  (63) 


The  rear  of  the  rod  is  undeformed  (du  =  d21  =0),  and  the  values  of  d12,  d22  are  given  in 
appendix  G.  Two  values  of  TBF  will  be  considered.  These  are  unbent  fins  and  the  rear  1  cal.  of 
fins  bent  uniformly  to  an  angle  of  0.02  radians. 

A  measure  of  the  flight  flexibility  of  a  missile  is  the  ratio  of  the  first  elastic  frequency  to  the  rigid 
missile  aerodynamic  frequency,  a  =  oox/coR  .  Calculations  of  <j>j  vs.  a  have  been  made  for  the  1-, 

3-,  5-,  and  7-element  rod.  The  results  for  all  but  the  1 -element  rod  are  practically  identical,  and 
the  curve  for  a  3 -element  rod  is  given  in  figure  4.  At  a  =  5 ,  the  aerodynamic  frequency  is  60% 
of  the  rigid  missile  aerodynamic  frequency;  and  even  at  ct  =  10  ,  it  is  10%  less  than  the  rigid 
missile  aerodynamic  frequency. 
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Figure  4.  <j>,  /  a>R  vs.  a  for  finned  missile. 

The  first  six  positive  frequencies  and  their  damping  rates  for  a  =  5  have  been  calculated  for  a 
1-,  3-,  5-,  and  7-element  rod.  These  results,  together  with  calculations  from  the  theory  of 
reference  (72),  are  tabulated  in  table  1.  We  identify  frequencies  which  differ  from  the  reference 
(72)  value  by  >10%  as  poor  values  and  mark  them  with  x’s.  Thus  the  1-element  rod  yields  one 
moderately  good  result  while  the  3 -element  and  5 -element  rods  yield  three  and  five  good  results, 
respectively.  For  this  elastic  problem,  n_j  —  1,  elastic  frequencies  are  well  calculated.  For  zero 
spin,  the  negative  frequencies  are  the  negatives  of  the  positive  frequencies. 

Because  the  3 -element  rod  predicts  the  lower  three  frequencies  quite  well,  all  of  the  calculations 
given  in  this  report  are  based  on  a  3-element  rod.  Most  of  these  have  been  repeated  for  a 
5 -element  rod,  and  almost  identical  results  have  been  obtained. 

In  figure  5,  f2  for  a  =  5,  TBF  =  0  is  plotted  vs.  <j)/coR  for  spin  near  the  first  elastic  frequency. 
Large  resonance  values  occur  at  <j)/coR  =  5.1  and  an  fj  line  for  pss  =6.0coR  is  also  shown.  This 
fj  line  has  three  intersections  with  the  f2  curve,  which  are  values  of  equilibrium  spin.  The 
intersection  near  pss/coR  is  denoted  as  pss/coR  ,  while  the  other  two  are  ple/coR  ,  p2e/coR  ,  and  all 
their  numerical  values  are  given  in  table  2. 
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Table  1.  Transient  frequencies  and  damping  rates. 


o  =  5  <j>  =  0 

Code 

k 

•k  W  R 

>-k/«R 

1  element 

1 

0.645 

-0.0530 

3  element 

1 

0.612 

-0.0528 

5  element 

1 

0.613 

-0.0528 

7  element 

1 

0.613 

-0.0528 

reference  (12) 

1 

0.614 

-0.0529 

1  element 

3 

6.146x 

-0.0855 

3  element 

3 

5.194 

-0.0724 

5  element 

3 

5.181 

-0.0718 

7  element 

3 

5.180 

-0.0717 

reference  (12) 

3 

5.184 

-0.0717 

1  element 

5 

20.31x 

-0.0644 

3  element 

5 

13.81 

-0.0518 

5  element 

5 

13.79 

-0.0523 

7  element 

5 

13.76 

-0.0518 

reference  (12) 

5 

13.74 

-0.0506 

1  element 

7 

— 

— 

3  element 

7 

30.00x 

-0.1018 

5  element 

7 

26.97 

-0.0914 

7  element 

7 

26.79 

-0.0881 

reference  (12) 

7 

26.69 

-0.0825 

1  element 

9 

— 

— 

3  element 

9 

53.37x 

-0.2060 

5  element 

9 

44.28 

-0.1389 

7  element 

9 

44.18 

-0.1462 

reference  (12) 

9 

43.84 

-0.1427 

1  element 

11 

— 

— 

3  element 

11 

100.44x 

-0.4708 

5  element 

11 

72.19x 

-0.1309 

7  element 

11 

65.88 

-0.1256 

reference  (12) 

11 

64.49 

-0.1202 

In  figure  6,  spin  histories  are  plotted  for  all  initial  conditions  equal  to  zero  except  <j>0  =  0,  5.4coR  . 
For  low  values  of  initial  spin,  lockin  occurs  at  the  first  elastic  frequency,  ple ,  and  large  strains 
result.  For  larger  values  of  initial  spin,  lockin  occurs  at  psse .  This  behavior  was  predicted  by  the 
simple  three-body  theory  of  reference  (11). 

It  is  very  unlikely  that  sufficient  fin  damage  occurs  to  produce  a  pss  in  excess  of  five  times  coR  . 
Because  it  is  much  more  likely  that  a  pss  near  coR  can  be  produced,  the  observed  inelastic 

deformation  is  probably  caused  by  lockin  at  the  aerodynamic  frequency,  and  this  case  is 
considered  in  detail  in  the  following  section. 
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Figure  5.  f,  vs.  <j)/ coR  for  T BF  =  0,  o  =  5,pss  =  6coR  . 


Table  2.  Equilibrium  spin  rates. 


0  =  5 

r 

*  BF 

0 

0 

0.02 

PSS/®R 

6.000 

0.800 

0.400 

PSSe/®R 

5.944 

0.788 

0.429 

Pl=/®R 

5.097 

0.590 

0.564 

Pss/®R 

5.393 

0.647 

0.627 

In  figure  7,  f2  is  plotted  for  spins  near  the  aerodynamic  frequency,  and  1',  is  shown  for 
pss  =  0.8cor  .  The  three  intersection  points  are  given  in  table  2.  Spin  histories  are  plotted  in 
figure  8  for  zero  initial  conditions  except  for  \j/0  =  0,  -  3  s'1 .  Here,  we  see  that  psse  spin-yaw 
lockin  occurs  for  zero  initial  angular  velocity,  and  ple  spin-yaw  lockin  for  an  initial  yaw  rate  of 
-3  s'1.  The  angular  motion  in  body-fixed  coordinates  is  plotted  in  figure  9.  For  lockin  at  psse , 
the  motion  is  smaller  than  0.09  rad,  while  for  lockin  at  ple ,  it  is  as  large  as  0. 15  rad.  The  motion 
of  the  forward  end  of  the  rod  in  body-fixed  coordinates  is  given  in  figure  10.  Motion  for  psse 
lockin  is  <0.17  (0.7  in),  and  after  1  s,  it  is  never  >0.09  (0.4  in).  For  ple  lockin,  it  exceeds  0.28 
(1.2  in). 
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Figure  6.  <j)/ cbr  vs.  time  for  Fbf  =0,  a  =  5,  pss  =  6.0coR,  <j)0/ coR  =  0,  5.4  . 
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0.05 
0 

-0.05 

-0.1 

Qb7z 

-0.15 

-0.2 

-0.25 

-0.3 

Figure  10.  Rod  forward  tip  motion  (qb7y  vs.  qb7z)  for  T BF  =  0,  a  =  5,  pss  =  0.8coR,  vj/0  =0,-3  s  '. 


-0.2  -0.15  -0.1  -0.05  0  0.05  0.1  0.15  0.2 


The  maximum  strain  on  the  rod  is  located  on  the  rod’s  surface  at  its  center.  In  terms  of 
3 -element  parameters,  the  maximum  strain  is 


s 


M 


(1/2) 


a:M°) 

dx2 


L;|l2q,  +  2LA 


(64) 


The  time  history  of  the  maximum  strain  is  plotted  in  figure  1 1  for  vj/0  =  0,  -  3  s'1 .  For  most 

metals,  yield  occurs  for  maximum  strains  in  excess  of  0.0015,  and  we  see  that  yield  does  occur 
for  resonance  lockin. 

The  two  regions  of  the  initial  angular  rate  plane,  i.e.,  the  \j/0  -0O  plane,  which  induces  either 

lockin,  can  be  detennined  by  a  number  of  trial  and  error  calculations,  and  the  boundary  curve 
between  these  regions  is  shown  in  figure  12.  Initial  conditions  outside  this  curve  will  cause 
resonance  lockin  while  conditions  inside  it  induce  psse  lockin. 

For  TBF  =  0  ,  the  first  resonance  peak  of  the  f2  curve  is  negative,  and  pss  must  be  greater  than 
the  aerodynamic  frequency  for  three  intersections  to  occur. 

As  is  shown  in  figure  13,  bent  fins  corresponding  to  TBF  =  0.02  - 1 1  <  x  <  -10  will  produce  a 
positive  peak  in  the  f2  curve  and  a  pss  less  than  the  aerodynamic  frequency  can  produce 
resonance  lockin.  The  plots  corresponding  to  figures  8-12  are  given  by  figures  13-19. 
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Figure  11.  Maximum  strain  ( em )  for rBF  =  0,  ct  =  5,  pss=0.8coR,  vj/0  =  0 ,  -3  s"1. 


9.  Summary 


The  application  of  FEM  to  the  calculation  of  the  motion  of  a  pitching,  yawing,  and  flexing 
missile  has  resulted  in  the  construction  of  a  family  of  linear  differential  equations. 

Computer  codes  have  been  developed  for  1,  3,  5,  and  7  elements  and  used  to  obtain  the  first  four 
elastic  frequencies.  The  3-element  code  has  given  excellent  values  of  the  first  two  elastic 
frequencies.  The  5-element  and  7-element  codes  provide  equally  good  values  for  the  first  four 
elastic  frequencies. 

Integration  of  these  differential  equations  has  demonstrated  resonant  spin  lockin  at  the 
aerodynamic  frequency  and  the  first  elastic  frequency.  Lockin  can  occur  near  the  design  steady 
state  spin  or  near  resonance  with  a  transient  frequency.  Calculations  have  identified  regions  of  a 
yaw  rate-pitch  rate  plane  that  determine  which  lockin  will  occur. 
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Figure  13.  fj  vs.  cj) / coR  for  rBF  =  0.02,  a  =  5,  pss  =  0.4g)r  . 
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Figure  15.  Angular  motion  (qbjy  vs.  qbiz)  for  Fbf  =0.02,  a  =  5,  pss  =  0.4coR,  0O  =0,  5  s'1. 
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Figure  16.  Rod  forward  tip  motion  (qb7y  vs.  qb7Z)  for  T BF  =  0.02,  a  =  5,  pss  =  0.4raR,  90  =  0  s'1. 
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Appendix  A.  Integrals 


A.2  Boundary  Terms 


xl 

x3 

I.= 

J  cfldx 

=  J  cfldx 

x0 

X2 

X1 

X3 

I3  = 

J  (x-Xj)cfldx 

I4  =  J(x-x2)  cfldx 

xo 

X2 

X1 

X3 

I5  = 

J  cf2dx 

l6  =  J  cf2dx 

x0 

x2 

X1 

X3 

I7  = 

J  (x_x1)cf2dx 

I8  =  J(X-X2)  Cf2dx 

x0 

X2 

X1 

X3 

I9  = 

J  (x-xj2  cfldx 

!,<>  =  J(x-x2)2cfldx 

x0 

X2 

1,1  = 

2I5-I3-XiIi 

Ij2  =  2I6  — 14  —  x2I2 

Il3  = 

2I7 -I9 -xxI3 

Il4  =  2Ig  — 110  —  X2I4 

X1 

x3 

IlD  = 

J  cDdx  +  CDbp 

I2d  —  J  cDdx 

x0 

x2 

X1 

x3 

I3D  = 

J  (x-Xj)cDdx 

I4D  =  J(x-x2)  cDdx 

xo 

X2 

xi 

X1 

llBF  — 

J  ^BF^BFdx 

I3BF  =  J(X_  Xl)  ^BF^BFdx 

xo 

xo 

xi 

xi 

-^5BF 

-  J  cf2FBFdx 

^7bf  =  J(x— xi)  cf2rBFdx 

x0 

xo 

(A-2) 
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A.3  Bent  Missile  Terms 


x2 

Jib  =  j*  cfi  ( +  rBF  )  dx 

x0 

x3 

^4B  =  1  [Cf2  FfiF  )  _  Cfl^EB  ]  xdx 

xo 

x3 

^2B  =  1  [Cf2  (^B  +  ^BF  )  _  Cf  l^EB  ]  dx 

xo 

x3 

^5B  =  J  CD  (^EB  —  dcB)dx  +  [^EB  (Xl)  ~  ^cB  ]  CDbp 

xo 

x3 

^3B  —  J  Cfl  (J^B  ^BF  )  xdx 

xo 

x2 

J6B  =(VL)  Jx5BPldx 

xl 

5cB=(VL)l5BPldX 

xi 

x2 

■^8B1  —  ~ad  J  rBP2dx 

xl 

x3 

IxB1=(l/L)J[pi8EB8EB+(p2adL)rBrB]  dx 

x3 

FB2  =  “6ad  J  [p2r BrB  ]  dx .  (A-3) 

*0  x0 


30 


Appendix  B.  Functions 


5c  =(1/L)J5EPldx 


^3 

J4(t)=  |(cf2r-cfI8E)  xdx 


Ji(t)  =  lcfirdx 


J5(t)  =  }cD(5E-Sc)  dx  +  [8E(x1)-8c]  cr 


J2  (t)  =  /  (cf2r-cfi5E)  dx  J6  (t)  =  (!/L)  J  x5EPidx 


J3  (t)  =  J  cfirxdx 


J8(t)  =  adj(r-2ici)r)p2dx. 


(B-l) 
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Intentionally  left  blank. 
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Appendix  C.  Generalized  Forces  and  Hermitian  Polynomials 


dF 

dQix  =  — JLdx  =  -g1cD(x)  dx 
dx 


dQiy  +idQlz  = 


dF  .  dF  ,  dF 

- +  1 +  (  -Ml  +  10 ) - 

dx  dx  ’  dx 


dx 


dQ2x  = 


dM„ 


dx 


dFy  .  dF 

— -  +  i — - 

dx  dx 


-6, 


dx 


dx 

dQ2y+idQ2z  =  |  x 

Qix=Fx 

Qiy  +  iQiz  =  Fy  +  iFz  +  (-V  +  i0)  Fx 
Q2x  =Mx 

Q!y+iQ,,=-i(My+iM,-iFA)-eM, 


dF  „  dM, 


-0- 


dx 


dx 


N,  =  l-3z2  +2z3 
N2  =  Lez(l-z)2 
N3  =  z2  (3-2z) 
N4  =  Lez2  (z-1) 


N(  =  6(z2  -z) 

N2  =Le(l-4z  +  3z2) 
=6(z-z2) 

N4  =  Le  (3z2  -2z). 


(C-l) 
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Appendix  D.  Connector  Coefficients  for  Equations  (39-41) 


D.l  Rod  Parameters 


^=Lejpi(x)  Nq(z)  dz 


:Le(zj+z) 


aj  =  LeJ  xpi  (x)  Nq(z)  dz 


Zj=(Xl/Le)  +  j-1 


^=Jp2(x)  Nq(z)  dz 


q  =  1,  2,  3,  4. 


(D-l) 


n.  +1 

For  central  element  r  =  — - 

2 


aic  -alr+5a3r  +  24Lea4r 


&2c  ~^2r  +Le^3r  +5^4, 


n.  +  3 

For  adjacent  element  r  =  — ^ — 


aL  -  5a[r  +  24Le'a2r  a2a  -  Lea|r  +  5a2r 


^3a  ^3r 


&4a  -  Kv 


aqc,  aqa,  bqc,  and  bqa  are  computed  from  aqj  and  bqj  in  same  manner  as  previously  shown. 


qj  qj 


ail  ai2  9 


n  =  4  L  =  L 


au  a  a14  a2c 

lc 


nt  =  8  Le  =  L/3 


an  an  a14  a21  a15  a31+alc+ala  a16  a41  +  a2c+a2a  an  a3a  ais  a4a 


nt  =  12  Le  =  L/5 


ai5  —  a31  +  ai2 


ai6  —  a41  a22 


ai7  “a32  +aic  +aL  aiS  -  a42  +a2c  +a2a  S19  “  aL  +  ai5 

/V  1  /V  1  /V  1  /v  1 

ai(10)  -a4a  +a25  ai(ll)  a35  ai(12)  a 45 
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a2n  andb2n  are  computed  from  a^,  a2qc,  a^a,  bf,  trc ,  and  Ir ;i  in  the  same  manner  as  previously 
shown. 
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Appendix  E.  Connector  Coefficients  for  Equations  (42^44) 


E.l  Rod  Terms 


fi=jcfl(x)  N;(z)  dz 
0 
1 

g!u  =  J[Cf2  (X)  NUZ)“Cfl(X)  Nq(Z)  Le]  dz 

0 

^  =  J[cfl(x)  Ni(z)] xdz 

o 

i 

g^i  =  j[Cf2(X)  NUZ)“Cfl(X)  Nq(Z)  Le]  XdZ’ 


?c  =  Le(zj+z) 

ZJ  =(Xl/Le)+ j-1 

q  =  1,  2,  3,  4 


(E-l) 


f , n ,  f2n,  gln ,  and  g2n  are  computed  from  fqj,  fq|,  gqj,  and  gq|  in  the  same  manner  as  in  appendix  D. 

E.2  Nonzero  Aerodynamic  Extension  Terms 


nt  =  4 


fai3=24L-e‘I2 
fa23  =24L-e1I18 


fal4=Il+5I2 

fa24  =  ^17  ^18 


gai3  =  -Ii  -  5I2  +  24L;1  (I6  - 14  )  gal4  =  I5  - 13  -  I2Le  +  5 (I6  - 14  ) 

ga23  =  — -^-17  —  ^I18  +  24Le  I16  ga24  =  I15  —  I18Le  +  5I16 


gal3  ^1 
ga23  _ ^17 


^a3  ~  ^ID 


2D  +  24Le  I4D 

da4  ^3DE 

nt  *4 

+  LeI2D  +51 

4D 

f 

xaln. 

=  I2 

f  =  I 

f 

=  1 

1a24  a17 

xa2n, 

a18 

ga14  =  ^5  _  ^3 

galf^-l)  ^ 2 

Sain, 

=  16- 

ga24  =  ^15 

gal(n,  — 1)  _ ^18 

Sain, 

=  Il6 

da4  =  ^3D 

^a(nt-l)  ^2D 

han, 

=  Cd- 

(E-2) 
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Appendix  F.  Differential  Equations  Coefficients 


F.l  Equation  44  m  =  1,  2 


Rm  =(md2/L)  aln 

R2„=(md2/L)  [a2n  +(adL)  b2n+xcaln] 

Su=(md)  V  +  (gld2/v)  c2 

S12  =  xcmd2 

Sin  =-(gld2/v)  (gln+galn) 

S2i=(gid2/V)  (c4-xcc2) 

S22  =It  +  md2x2 

S2n  =-(gid2/Y)  (g2n+ga2n)-xcSln 

Tn  =(gid)  (c!-Cd) 

T,2=(md)  V 

Tln=-(g1d)  (fln+faln) 

T2i=(g,d)  (cj-x^) 

T2n=-(gld)  (f2n+fa2n-CDair 

s2„  =-2(  md2)  adb2n 

A 

II 

1 

HH 

X 

tDi=(<i)2-4)  scB 

b)2  (46B  Al) 

*-ai  =4ib  (  ‘M/v)  J2B 

^  A2  =  A  +  *  (  ‘M/v)  J4B 

F.2  Equation  44  m  =  3,  4. .  .nt 

R™,=HJ/L)km 
S„,=(md)  Vd.,- 

(gid2/v)  (gml+gaml)  Sm2  =  R2m 

Sm„=-(gid2/v)  (gmn+gainn)+  (md2/L)  ( 2co2 /coj )  kcmn 

Tml  =  -fed)  (fml  +  faml  +  CDalmL-1 )  Tm2  =  mdVa^L-1 

Trrai=-(g,d)  (fm„+fanm)  +  (md2/L)  ®ocmn  +  42md2adbmn 

SL„  =  -4(  md2)  T*2  =  2md2adb2m 

t,L  — (md2/L)  (2co2/co1)  kCmn 

*mD  =  -  4)  (kBm  /L  +  2a  Am  )  "  Am 

=  ^B3  ^1B  i(  ^d/V^)  ( gB3  "^5B  ) 

A  =  ^B4  ^3B  i(  ^d/v^)  (gB4  ^7B  ) 
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tmA  =fBm+i(  ^/V)  gBm  4  <  111  <  II,  -  1 

'  t(nt-l)A  =  +  ^B  (X2  )  +  *  (<j)d/V )  (g^)  +  1^ B  (x2  )) 

tn.A  =fBnt  +  ^B  (X2  )  + 1  («j»d/V)  (§Bnt  +  I8rB(X2)) 

For  1^=1  nt=4  m  =  3,  4 
tmD  =(<l>2-4)  kBm  /L 

t3A  =  L  +  IlB  +  5I2F B (X2 ) +  24Le%r B  (  X2  ) 

+  i(4»d/v)  ( gBic  +  I5B  +  5I6rB  (x2 )  +  24L-e1I8rB  (x2 )) 

^4A  =  ^B2c  ^3B  +Lel2rB  (X2)  +  5I4FB  (X2) 

+  i(cj)d/v)  (gB2c  +  I7B  +  LeI6r B  (x2 )  +  5Isr B  (x2 )).  (F-2) 


dx 


F.3  Quadratic  Spin  Equation  (61) 

J7  =  -(i/L)  J  Lpi8eSe  -  p2adL(r  +  4i<j)r  +  iq2)  F 

X1 

nt  nt 

Qd  =  -(i/L)  Z Z [(amn  -  adLbmn )  qmqn  +  8i(j)adLbmnqmqn ] 

m=3  n=3 

-  (i/L)  e_i+  Z  [(aBm  -  adLbBm )  [qm  +  4>2qm]  +  8icj)adLbBm  (qn  -  i(j)qm ) 

m=3  L 

Jx  =Ix+IxB-(md2/L)  Re(4adLZ^bmnqmqn+e^Z(aBm-3(adL)bBm)  qmj 


m=3  n=3 


m=3 


^xB  md  (IXB1  IXB2  )• 


(F-3) 


40 


Appendix  G.  Flexing  Motion  Finite  Element  Method  (FEM)  Quantities 


G.l  Kinetic  and  Potential  Energy  Coefficients 


apqj=Lejpi(x)NpNqdz 


,  J.  ,  ,  dNn  dNn 

b»=L'/p=(x)^7l7dz 


£  -L-3fp  (x)^X£Xdz 
pcu  eJp3^  dz2  dz2 


(G-l) 


n  •  + 1 

2x2  central  element  matrix  r  =  — - 

2 

anc  =  3-iir  +  5 (a13r  +  a31r )  +  24Le'  (a14r  +  a41r )  +  25a33r  +  1 20LC'  (a34r  +  a43r )  +  (24)  a44rLe 
ai2C  =  ai2r  +  Lea13r  +  5  (a14r  +  a32r )  +  24Le‘a42r  +  5Lea33r  +  25a34r  +  24a43r  +  120Le'a44r 
a2ic  =  a2ir  +  5  (a23r  +  a41r )  +  24Le1a24r  +  Lea31r  +  5Lea33r  +  24a34r  +  25a43r  +  120Le1a44r 
a22c  =  a22r  +  Le  (a23r  +  a32r )  +  5  (a24r  +  a42r )  +  L;a33r  +  5Le  (a34r  +  a43r )  +  25a44r 

n  •  +  3 

4x4  adjacent  element  matrix  r  =  — 


ana  -  25aiir  +  120Le‘  (a12r  +  a21r )  +  (24)  a22rLe 
ai2a  —  5Leallr  +  24a21r  +  25a12r  +  120Le  a22r 


a2ia  5Leallr  +  25a21r  +  24a12r  +  120Le  a22r 
^22a  =  Leallr  +  5Le  (a12r  +  a21r )  +  25a22r 

^13a  =  ^^13r  ^4a23rLe  3 

a[4a  =5a14r+24a24rL;1  a 

a23a  =  Lea13r +5a23r  a 

^24a  =Le^l4r  +  5^24r  a 


^31a  ^^31r  24a32rLe 

^41a  —  ^^41r  +  24a42rLe 
a32a  =  Lea31r  +  5a32r 
a42a  =  Lea41r  +  5a42r 


a33a  a33r 


a43a  a43r 


a34a  a34r 


a44a  a44r 


bpqc,  b  ,  Cpqc ,  and  cpqa  are  computed  from  bpqj  and  cpqj  in  the  same  manner  as  previously  shown 


pqi  pqj 
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n,=4 


nt  =  8 

Le  “  L/3 


a33 

=  aill 

a34 

=  ai21 

a35 

—  ai31 

a36 

II 

£ 

a43 

—  a211 

a44 

—  a  22 1 

a45 

—  a231 

a46 

—  a241 

a53 

—  a311 

a54 

=  a321 

a55 

=  a331  +  allc  +  alla 

a56 

—  a341  +  a12c  +  a12a 

a57 

ai3a 

a58 

ai4a 

a63 

—  a411 

a64 

—  a421 

a65 

—  a431  +  a21c  +  a21a 

a66 

—  a441  +  a22c  +  a22a 

a67 

a23a 

a68 

a24a 

a75 

—  a31a 

a76 

—  a32a 

a77 

—  a33a 

a78 

—  a34a 

a85 

=  a41a 

a86 

“  a42a 

a87 

=  a43a 

a88 

=  a44a  other  amn 

's  are  zero 

n,  =12 
L,  =  L/5 


a33  “  aill 

a34  =  ai21 

a35 

“  a  13 1 

a36  = 

ai41 

P 

II 

N?> 

a44  =  a  22 1 

a45 

=  a231 

a46 

5  24i 

a53  —  a  3 1 1 

a54  —  a321 

a55 

—  a331  +  a112 

1 

a341  +  a  1 22 

a57  — 

ai32 

a58  —  a  142 

a63  =  a411 

a64  =  a  42 1 

a65 

=  a431  +  a212 

a66  = 

a  44i  +  a222 

a67  = 

a232 

a68  =  a242 

a75  —  a312 

a76  —  a322 

a.  77 

=  a332  +  allc  +alla 

a78  — 

a342  ai2c  +  ai2a 

a  79  — 

ai3a 

a 7(10)  _  ai4a 

a85  =  a412 

a86  =  a422 

a87 

=  a432  +a21c  +  a21a 

a88  = 

a442  +  a22c  +  a22a 

a89  = 

a23a 

a8(10)  =  a24a 

a97  =  a31a 

a98  =  a32a 

a99 

=  a33a  +ail5 

a9(10)  _ 

:  a34a  +  a  125 

a9(ll)  _ 

=  ai35 

a9(12)  —  a  145 

(10)7  _  a41a 

a (1 0)8  _  a42a 

a(10)9 

—  a43a  +  a215 

a(io)io  “ 

'  a44a  a225 

a(io)ii  “ 

“  a235 

a<  10)12  —  a245 

L(1 1)9  =  a  3 1 5 

a(ll)10  =  a325 

a(ll)ll 

=  a335 

a(l  1)12  _ 

:  a345 

(12)9  =  a415 

a(l2)10  —  a425 

a(l2)ll 

=  a435 

a(12)12 

=  a445  other  a 

^  's  are  zero 
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nt  =  16 
L,  =  L/7 


a33 

—  aill 

a34  = 

ai21 

a35  ai31 

a36  : 

“  ai41 

a  43 

=  a211 

a44  — 

a221 

P 

Ji¬ 

ll 

a46 

=  a241 

a53 

=  a311 

a54  ~ 

a321 

a55  ~  a331  ai  12 

a56 

—  a341  +  a122 

a57  ' 

“  ai32 

a58  ~  ai42 

a63 

—  a411 

a64  — 

a421 

a65  —  a 43 1  a 2 12 

a66 

—  a441  +  a222 

a67  ' 

—  a232 

a68  —  a242 

a75 

=  a312 

a76  = 

a322 

a77  —  a332  ai  13 

a78  : 

—  a342  ai23 

a79  = 

=  ai33 

a7(10)  —  a  143 

a85 

—  a412 

aS6  = 

a422 

a87  —  a432  a213 

a88  : 

—  &442  +  a223 

a89  = 

=  a233 

a8(10)  —  a243 

a97 

=  a313 

a98  = 

a323 

a99  =  a333  +ailc  +aila 

a9(10) 

a343  +  ai2c  +  ai2a 

a9(ll)  ' 

=  ai3a 

a9(12)  —  ai4a 

a ( 10)7 

—  a413 

a(10)8  “ 

'  a423 

a(10)9  —  a433  a21c  a21a 

a(io)io 

—  a443  +  a22c 

a22a 

a(io)n 

—  a23a 

a(10)12  —  a24a 

a(ll)9 

=  a31a 

a(n)io  = 

'  a32a 

a(l  1)1 1  =  a33a  +ail6 

a(l  1)12 

—  a34a  ai26 

a(ll)13 

II 

OS 

a(l  1)14  —  ai46 

a(12)9 

—  a41a 

a(  12)10  “ 

'  a42a 

a(12)ll  —  a43a  a 2 1 6 

a(12)12 

a44a  a226 

a(12)13 

a236 

a(12)14  —  a246 

a(  13)1 1 

=  a316 

a  ( 13  )12  ““ 

'  a326 

a(13)13  =  a336  +ail7 

a(13)14 

—  a346  ai27 

a(13)15 

—  ai37 

a(13)16  =  ai47 

a(14)ll 

=  a416 

a  ( 14)12  = 

"  a426 

a(l4)13  —  a436  a 2 17 

a(14)l4 

—  a446  ^227 

a(14)15 

=  a237 

a(14)16  ~~  a247 

a(l5)13 

=  a317 

a(15)14  = 

=  a327 

a(15)15  —  a337 

a(  15)16 

—  a347 

a  ( 16)13 

=  a417 

a<  16)14  = 

=  a427 

a(16)l5  —  a437 

a(16)16 

—  a447 

other 

amn  's  are  zero 

b  ,  and  Cllm  computed  from  b^,  bpqc,  bpqa,  c^,  Cpqc,  and  c^a  in  the  same  manner  as  previously 


shown. 

kmn  =amn+(adL) 

G.2  Kinetic  Energy  FEM  Vectors 

1 

aBpj  =  Lejpi  (x)  8eb(x)  Np(z)  dz  x  =  Le(Zj+z) 

0 

bBPj=Jp2(x)  Fb(x)  Np(z)  dz  Zj  =(xj/Le)  + j-1.  (G-2) 

0 

aBm  and  b.Bm  are  computed  from  aBpj  and  bBpj  in  the  same  manner  as  aln  was  computed  from  af, 

kBm  =  aBm  +  (adF)  bBm  • 
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x  =  Le  (Zj  +z) 


G.3  Rod  Aerodynamic  Force  Terms 

fpqj=jcfl(x)  NqNp(z)  dz 
0 

gpqj=  j[cf2(x)  Ni(z)-cfl(x)  Nq(z)  Le]  Np (z)  dz  Zj  =(Xl/Le)  + j-1 
0 

P,q  =  l,  2,  3,  4.  (G-3) 

fmn  and  grnn  are  computed  from  fpqj  and  gpqj  in  the  same  manner  as  arnn  was  computed  from  apqj. 


fpji=-Lcjcfi(x)  Np  (z)  dz 

0 

x  =  Le(zj  +  z) 

1 

gPJl=-LeJCf3(x)  Np  (z)  dz 

o 

zJ=(xi/Le)  +  j-1 

(x)  =  2cf2  —  xcfl 

p  =  1,  2,  3,  4 

fBpj  =LeJcfi(x)  [rB(x)  +  rBF(x)]  Np  (z)  dz 
0 
1 

gBPj=Lejcf2(x)  [rB(x)  +  rBF(x)]  (x)  Np (z)  dz 
0 

fmi  ’  g m i  ’  f'linP  and  gBm  are  computed  from  fpjl,  gpjl,  fBpj,  and  gBpj  in  the  same  manner  as  aln  was 
computed  from  a f.  All  other  fml,  gml,  fBm,  and  gBm’s  are  zero. 

G.4  Nonzero  Aerodynamic  Extension  Terms 

n,=4 

f.j,  =-[!,  +5I2  +  24L;[I4]  f,M  =  24L;[  [5(I2  +  I,d)+  24L;‘I4] 

f,«  =  [I,  +  251.  +  241.  „  +  120L;1  (I4  +  I4D )]  f,41  =  -  [I,  +  L,I2  +  5I„] 

f,«  =  [cDlp  +  241.  +  25I.d  +  120L;1  (I4  +  I4D)]  f.„  =  [i,  +  5L,  (I,  +  I2D)  +  25(1,  +  I4D)] 

f.Rj  = 1  l'i,('x,)+ si.r B (x, )  +  24L;‘i4r B (x, ) + iIBF  f„4  =  isrB (x, )  +  5i„rB (x, ) + 24l  'ij  b  +  i3BF 
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nt  *  4 


f.3i  =  -1: 
f  =  -I 


fa34  -  I, 
f  =T 

±a(nt-i)nt  2 


f*  —  c 

1a43  ^Dbp 

f  =  T 

ant(nt-l)  2D 


f.41  =  -J3 

f  =-T 

a(n,)l  A4 


f  =1 

a  44  a3 

fan.n,  =  ^4  ^4D 


faB3  =IirB(Xl)  +  IlBF  faB4=I3rB(Xl)  +  I3BF  faB(„,-l)  =  ^B  (X2  )  faBn,  =  U^B  (  X2  ) 

ero 

n  =4 


other  fAnm  's  are  zero 


ga31--[ln+5I12  +  24L;I14]  ga33=-Ll1+25I2  +  120L;1(2I4-I6)  +  (24)2L-e2(l10-I8) 

ga34  =  ■ -  [I3  -  Is  +  5LeI2  +  24  (2I4  - 16 )  +  25  (l10  - 18 )] 

ga4i  =-[li3+  LcI12  +  5I14 ]  ga43  =  - [l3  +  5LeI2  +  24 (2I4  - 16 )  +  120L;1  (l10  - 18 )] 

ga44  =  -[h>  “  h  +  Kh  +  5LC  (2I4  - 16 )  +  25  (l10  - 18 )] 
gaB3  =  ^B  (X1  )  +  5I6r B  (X2  )  +  241^1^ B  (x,  )  +  I5BF 
gaB4  —  ^7^B  (xi ) + Lei6rB  (x2 ) + 5i8rB  (x2)+ i7BF 

nt  *  4 


ga31  bl 
ga41  —  _  ^13 

=  -I, 


ga33  b 

ga43  =  — ^3 
a(n„-l)(nn-l)  ^2 


5a34  b  ^3 
ga44  =  1 7  _  ^9 


^a(n„-l)(nn)  "6  "4 


=  W 

8a(n,)(nt)  _  ^  ^ 


10 


®a(nn-l)l 

§a(nt)l  —  ^14  §a(nt)(n,-l)  —  ^4 

gaB3  =  ^5^8  (Xl)  +  IsBF  gaB4  =  (  X1  )  +  ^7BF  gaB(n,  -l)  =  (  X2  )  gaBn(  =  ^8^8  (X2  ) 

other  ga 


s  are  zero  . 


(G-4) 
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Intentionally  left  blank. 
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Appendix  H.  Finned  Missile  Parameters 


Pi  P2  P3  1 


=  0 


L  =  20 
d  -  0.35  ft. 
m  =  3.50  slug 
Ix=  0.054  slug-ft2 

It=  14.318  slug-ft2 
Cfi  =4(1 1-x) 

=  e2(*-10) 

=  -(2/7t)(l5  +  3x)aL 
=  (12A)a  L 
cf2  =  2(1 1-x)2 

=  2  +  0.5(l-e2(x~10)) 

=  2.5 +  (l/37i)  (15  +  3x)2  aL 
=  2.5  —  (12/ti:)  (6  +  x)  aL 
cD  =(0.30)  (11-x) 

=  0 

CDbp=0.14 

di2  =10~3 
CNa=9-7 
CMq=-980 


V  =  6000  ft/s 
p=  0.002  slugs/ft3 


X01  X23  1 

as=  1100  ft/s 


10  <  x  <  1 1 

-5  <  x  <  10 
-7  <  x  <  -5 
-1 1  <  x  <  -7 


10  <  x  <  1 1 

-5  <  x  <  10 
-7  <  x  <  -5 
-1 1  <  x  <  -7 


10  <  x  <  1 1 
x  <  10 


c* 

d22 

c 

'-'Ma 

r 

^Md 


=  -18 

=  -0.25x10  5 
=  -34.4 
=  -190 
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Intentionally  left  blank. 
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List  of  Symbols,  Abbreviations,  and  Acronyms 

cfj  (x)  aerodynamic  force  distributions  functions 

d  rod  diameter 

E  Young’s  modulus 

F  Fy  +  iFz  complex  transverse  aerodynamic  force 

a„  (16L)”1 

gi  g,  pV=S/2 

I  (d)4  ff  y2dydz-(d)4f[z’dydz  ,  area  moment  of  rod 

Ix  axial  moment  of  inertia  of  projectile 

It  transverse  moment  of  inertia  of  projectile 

k  beam  damping  coefficient 

F  rod  length  /  rod  diameter 

m  projectile  mass 

n  ■  number  of  rod  elements 

nt  2nj  +  2 

p  projectile  spin 

q,  []  +  ia ,  complex  angle  of  attack  of  central  disk  (nsc) 

qle  \|/  -  i0  complex  yaw  and  pitch  of  central  disk  (nsc) 

q2  \j/  -  i0  complex  yaw  and  pitch  rate  of  central  disk  (nsc) 

qn  n  =  3 .4 . . .  nt  FEM  connectors  (nsc) 

qbn  n  =  3 .4 . . .  nt  FEM  connectors  (bfc) 

Fe  L/nj 

S  7rd2/4 
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T 


total  kinetic  energy 

kinetic  energy  of  disk  at  x  with  thickness  dx 
magnitude  of  projectile  velocity 


Td 

V 

Xj ,  x2  location  of  beam  ends 

x01,x23  dimensionless  length  of  fore  and  aft  aerodynamic  extensions 
xc  axial  location  of  center  of  mass 

a  angle  of  attack  of  central  disk 

P  angle  of  sideslip  of  central  disk 

T  — ^ ,  complex  cant  of  disk 

dx 

sM  maximum  strain  of  rod 

X, 

Sc  If1  J  SEdx ,  lateral  location  of  missile’s  cm 

xi 

5e  8e  +  i§Ez ,  lateral  displacement  of  disk  (nsc) 

(j)  roll  angle 

cj)k  frequency  of  k-th  mode 

kk  damping  of  k-th  mode 

p  air  density 

pj  axial  variation  of  mass,  moment  of  inertia,  elasticity 

a  cOi/c^R 

®o=  E0I0(L/md) 

C0j  lowest  elastic  frequency  of  beam  in  vacuum 

coR  rigid  projectile  zero-spin  frequency 

F  =  ( Fx ,  Fy ,  F7  j ,  aerodynamic  force  exerted  on  missile  (nsc) 

M  =  (Mx,My,Mz),  aerodynamic  moment  exerted  on  missile  (nsc) 
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Re{z}  real  part  of  z 
Tm  { zj  imaginary  part  of  z 

Carat  superscript  denotes  quantity  for  a  single  element 
Tilde  superscript  denotes  elastic  parameter  for  bent  missile 
B  subscript  denotes  parameter  for  bent  projectile 
BF  subscript  denotes  bent  fin  parameter 
E  subscript  denotes  an  elastic  coordinate  parameter  (nsc) 

T  subscript  denotes  trim  motion  parameter 
b  subscript  denotes  an  elastic  coordinate  parameter  (bfc) 
(bfc)  body-fixed  coordinates 
(efc)  earth-fixed  coordinates 
(nsf)  non-spinning  coordinates 
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07806-5000 


4  COMMANDER 
US  ARMY  ARDEC 
SMCAR  FSE 
A  GRAF 
D  LADD 

E  ANDRICOPOULIS 
K  CHEUNG 

PICATINNY  ARSENAL  NJ 
07806-5000 

6  COMMANDER 
US  ARMY  ARDEC 
SMCAR  CCL  D 
F  PUZYCKI 
D  CONWAY 
D  DAVIS 
K  HAYES 
M  PINCAY 
W  SCHUFF 

PICATINNY  ARSENAL  NJ 
07806-5000 

3  DIRECTOR 

US  ARMY  RESEARCH  OFC 
G  ANDERSON 
K  CLARK 
T  DOLIGOWSKI 
PO  BOX  12211 

RESEARCH  TRIANGLE  PARK  NC 
27709-2211 

2  COMMANDER 

US  NAVAL  SURFC  WEAPONS  CTR 
CODE  DK20  MOORE 
CODE  DK20  DEVAN 
DAHLGREN  VA  22448-5000 

2  COMMANDER 

WHITE  OAK  LABORATORY 
US  NSWC  APPLIED  MATH  BR 
CODE  R44  PRIOLO 
CODE  R44  WARDLAW 
SILVER  SPRING  MD  20903-5000 

1  COMMANDER 

US  ARMY  AVN  AND  MIS  CMND 
AMSAM  RD  SS  AT 
W  WALKER 

REDSTONE  ARSENAL  AL 
35898-5010 
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COMMANDER 

US  AIR  FORCE  ARMAMENT  LAB 
AFATL  FXA 
B  SIMPSON 
G  ASATE 
R  ABELGREN 
G  WINCHENBACK 
EGLIN  AFB  FL  32542-5434 

DIRECTOR 

SANDIA  NATIONAL  LAB 
W  OBERKAMPF 
W  WOLFE 
DIVISION  5800 
ALBUQUERQUE  NM  87185 

DIRECTOR 

LOS  ALAMOS  NATIONAL  LAB 
MS  G770  W  HOGAN 
LOS  ALAMOS  NM  87545 

DIRECTOR 

NASA  AMES  RESEARCH  CTR 
MS  258  1  L  SCHIFF 
MOFFETT  FIELD  CA  94035 

DIRECTOR 

NASA  LANGLEY  RESEARCH  CTR 
M  HEMSCH 
LANGLEY  STATION 
HAMPTON  VA  23665 

MASSACHUSETTS  INST  OF  TECH 
DEPT  OF  AERONAUTICS  AND 
ASTRONAUTICS 
E  COVERT 

77  MASSACHUSETTS  AVE 
CAMBRIDGE  MA  02139 

ARROW  TECHLGY  ASSOC  INC 
R  WHYTE 
PO  BOX  4218 

BURLINGTON  VT  05401-4218 

C  H  MURPHY 
PO  BOX  269 

UPPER  FALLS  MD  21156 

W  H  MERMAGEN 
4149  U  WAY 

HAVRE  DE  GRACE  MD  21078 


1  DR  W  STUREK 

3500  CARSINWOOD  DR 
ABERDEEN  MD  21001-1412 

1  OREGON  STATE  UNIV 
DEPT  OF  MECH  ENGR 
DR  M  COSTELLO 
CORVALLIS  OR  97331 


ABERDEEN  PROVING  GROUND 

4  COMMANDER 
US  ARMY  ARDEC 
SMCAR  DSD  T 
R  LIESKE 
F  MIRABELLE 
J  WHITESIDE 
J  MATTS 
APG  MD  21005 

26  DIR  USARL 

AMSRD  ARL  Cl  H 
C  NIETUBICZ 
AMSRD  ARL  HR  SC 
D  SAVICK 
AMSRD  ARL  SL  BE 
A  MIKHAIL 
AMSRD  ARL  WM 
E  SCHMIDT 
AMSRD  ARL  WM  B 
A  HORST 

AMSRD  ARL  WM  BA 
W  D  AMICO 
B  DAVIS 
T  HAWKINS 
AMSRD  ARL  WM  BC 
M BUNDY 
G  COOPER 
W  DRYSDALE 
J  GARNER 
B  GUIDOS 
BOSKAY 
P  PLOSTINS 
JSAHU 

K  SOENCKSEN 
P  WEINACHT 
S  WILKERSON 
AMSRD  ARL  WM  BD 
T  MINOR 
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AMSRD  ARL  WM  BF 
HEDGE 

AMSRD  ARL  WM  T 
B  BURNS 

AMSRD  ARL  WM  TC 
R  COATES 
WDEROSSET 
RMUDD 

AMSRD  ARL  WM  TD 
E  RAPACKI 
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1  DEFENSE  RESEARCH 

ESTABLISHMENT  VALCARTIER 
DELIVERY  SYSTEM  DIVISION 
A  D  DUPUIS 
2459  PIE  XI  NORD 
VAL  BELAIR  QUEBEC 
CANADA 


56 


